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We present a model for a biological reactor in which bacteria tend to aggregate 
in floes, as encountered in wastewater treatment plants. The influence of this 
flocculation on the growth dynamics of the bacteria is studied. We argue that 
a description in terms of a specific growth rate is possible when the flocculation 
dynamics is much faster than the other processes in the system. An analytical 
computation shows that in this case, the growth rate is density-dependent, i.e., 
depends both on the substrate and the biomass density. When the flocculation 
time scale overlaps with the other time scales present in the system, the notion 
of specific growth rate becomes problematic. However, we show numerically 
that a density-dependent growth rate can still accurately describe the system 
response to certain perturbations. 



1 Introduction 



Biological reactors are commonly used to remove pollutants from wastewaters. One stan- 
dard technology is the two-step Activated Sludge Process (ASP). Both in the reaction 
and the settling tank, bacteria naturally aggregate and form floes. It is well known - but 
poorly understood - that the flocculation is strongly dependent on the mass and volu- 
metric loading rates: the higher the loading rate, the higher the risk of decreasing sludge 
floe formation and settling capacity. In order to optimize this bioprocess, it is therefore 
important to better understand the flocculation phenomenon. 

Mathematical modelling has turned out to be a valuable tool in the study of Wastewater 
Treatment Plants (WWTP). To coordinate efforts in the development of efficient design 
and control tools, the Activated Sludge Model No. 1 (ASM1) was proposed in the late 
eighties. This model describes the different biological processes (e.g., chemical oxygen 
demand (COD) removal, (de)nitrification and phosporous removal) in detail. Its core 
consists of the mass balance equations, including the reaction kinetics as a function of the 
limiting substrates, which read in their simplest form, 

-7- = h(s)x — Dx 
At K ' 

As 

— = -h(s)x + D(s iu -s), (1) 

where x is the biomass concentration, s the substrate concentration, h(s) the specific 
growth rate, D the dilution rate and s ln the substrate concentration in the inflow. 
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Although a number of improvements have been investigated for the reaction dynamics (e.g., 
ASM2, ASM2d and ASM3), the model part for the floe formation and settling remains 
the weakest part. This modelling problem (see ^1 for reviews) has been studied from 
different perspectives. Population Balance Models (PBM) describe the floe aggregation 
and breakage and allow to compute the floe size distribution as a function of time 01 
E|. Computational Fluid Dynamics (CFD) simulators describe the hydrodynamics in the 
clarification tank and try to predict the settling properties of the floes [B]. Individual- 
Based Models (IBM) take both physico-chemical and biological processes into account at 
the level of a single floe 0. 

These modelling approaches have in common a high-dimensional parameter space. Al- 
though these parameters can be identified from experiments, the resulting model is often 
too complex to provide insight in the governing mechanisms. Moreover, to compute, for 
instance, the settling properties of the ensemble of interacting floes in the clarifier, one has 
to combine a CFD with a PBM approach, which leads to even more intricate models. 

Instead of starting from performant simulators, we propose to take the simple model 
as a point of reference. In particular, we investigate how these equations are modified 
when the biomass is organized in floes. To address this question, we propose a PBM-likc 
model where both the floe interactions (as in standard PBM) and the bacteria growth are 
included. This qualitative model is sufficiently transparent to be manipulated analytically. 
Our approach is primarily intended to model the floe dynamics in the reaction tank, where 
both physico-chemical and biological processes have to be taken into account. Nevertheless, 
our model can also be useful to check the common assumption of PBM that biological 
growth can be neglected in the settling tank. 

Our analysis naturally leads to an effective model of the form, 

h(s, x)x — Dx 

-h[s,x)x + D(s ia -s). (2) 

Note that the specific growth rate h(s, x) depends both on the substrate concentration s 
and the biomass concentration x, in contrast with the substrate-dependent growth rate 
h(s) of model 0J. The specific growth rate h(s,x) is called density-dependent. In fact, 
based on the original work by Arditi and Ginzburg 8 , density-dependent growth rates were 
recently proposed to describe bioreactor kinetics more accurately From an ecological 
point of view, this change has important consequences, as it allows microorganisms to 
coexist in a medium where classical, i.e. substrate-dependent, models predict extinction 
by wash-out. 

It should be noted that this work is not the first to study the influence of a heterogeneous 
biomass structure on the growth rate (see, for example, |10l lllj). However, we present 
here, to the best of our knowledge, an original derivation of an effective model with density- 
dependent growth dynamics, starting from a PBM description including bacterial growth. 

This paper is organized as follows. In Section 2, we introduce the bioreactor model. The 
different phenomena, including bacterial growth, floe aggregation and breakage, and hy- 
drodynamics, are discussed. In Section 3, we perform an analytical study of the model, 
under the hypothesis that the time scale associated with the floe interactions is much 
shorter than the other processes present in the system. We show analytically how this hy- 
pothesis leads to a density-dependent growth rate. In Section 4, we present some numerical 
computations, that go beyond the hypothesis of separate time scales. 



2 Modelling flocculation of growing bacteria 

We start by introducing some notation. Consider a bioreactor in which a biomass grows 
on a substrate. The density of the biomass is denoted by x, the density of the substrate by 



da; 
df 
ds 
dt 
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s. The biomass consists of bacteria which naturally aggregate in floes. A floe containing n 
bacteria will be denoted by F n . Define u n as the density of floes of size n. Expressing the 
densities x resp. u n as the number of particles (bacteria resp. floes) per unit of volume, 
we have 



nu n . (3) 



We assume the reactor to be perfectly mixed. As a consequence, all floes have the same 
access to the substrate. However, we will take into account that the bacteria at the surface 
of the floes have easier access to the substrate than the bacteria inside the floes. 

The dynamics of the floe densities u n is given by 

du n ( du n \ ( du n \ 

~dT = \~dTI ~ n \~dTi ' ^' 

V / bacterial growth V / floe interaction 

The second term in the right-hand side represents the bacteria disappearing in the effluent 
of the reactor with dilution rate D. The other two terms, due to the bacterial growth and 
the floe interaction, are now described in more detail. 



Dynamics of bacterial growth 

The only bacterial growth present in our model is through cell division. As a bacterium 
present in a floe of size n divides, we assume the daughter bacteria to stick to the floe, 
which will then consists of n + f bacteria. This growth can be written as 

F n — ■> F n+ \ with reaction rate h n (s). 

Note that the growth rate of a floe depends on its size n and on the substrate density s. 
To describe that the substrate has immediate access to the outer shell of the floe, whereas 
the inner bacteria can be deprived from the substrate, we assume the dependency on n 
and s to be 

h n {s) = h(s)n a , (5) 

with the exponent a < 1. When the access to the substrate is not limited, a = 1; when 
it is, a < 1 and its value corresponds to the surface-to-volume ratio of the floes. For 
spherical floes, a = 3. As floes are known to have some fractal form, also other exponents 
are possible. 

As will become clear from the computations, our results do not critically depend on the 
function n 1— ► h n (s). A more complicated sublinear dependence than J^} can easily be 
handled, taking the hydrodynamics around the floe into account ^2], or using the results 
of an individual-based model 7 . As our treatment is qualitative, and the precise behaviour 
of h n (s) therefore of secondary importance, we will restrict our attention to floe growth 
functions of the form (j^J . 

The dynamics corresponding to bacterial growth is 

-hi(s)ux 

ft n _i(s)M„_i - h n (s)u n , n>2. (6) 

Indeed, a growth event F n — > F n+i corresponds to the consumption of a floe of size n and 
the production of a floe of size n + 1. Mass action kinetics are assumed for this reaction. 



dui\ 
dt J 

du r 
~~dl 



bacterial growth 



bacterial growth 
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Aggregation-breakage dynamics 



The floe interactions we consider are the aggregation of two floes to form one bigger floe 
and the breakage of one floe into two smaller ones. As equations (@J are continuous in 
time, processes involving three or more floes are implicitly included. The floe interactions 
can be written as 

F m + F n — ► F m+n with reaction rate a m ,ni 

Fm+n — * F m + F n with reaction rate b m ^ n . (7) 

The reaction rates are symmetric in their arguments, i.e., a m , n = a n<m and 6 mi „ = b n>m . 
Many studies have been carried out to derive theoretical relations between these coefficients 
[131 or to identify them from experiments 0E1E]- We note already that our analysis, 
in the first place qualitative, will not need explicit expressions for the reaction rates a TO) „ 
and b m . n . 

The part of the dynamics corresponding to the floe interactions is 

L"f J oo 



dt 



— ^ ®m,n— m'U'm'U'n— m ^ (1 ~t~ &rn,n)Q j rn,n f Urn'U>n 
floe interaction m—l m—l 

oo J 

~t~ ^ ^ (1 ~t~ ^m.n^)bm.n^m+n ^ ^ b mn —mV"n-> (8) 
■m—l m — l 

where [x\ denotes the largest integer smaller than x, and 

J 1 if m = n, 
I otherwise. 

The first term corresponds to the aggregation of two floes to form a floe F n . The second 
term corresponds to the aggregation of a floe F n with another floe. The third term 
corresponds to the breakage of a floe into two floes, one of which has size n. The fourth 
term corresponds to the breakage of a floe F n into two smaller ones. 

The equations JHJ are identical to the so-called population balance models (PBM), first 
introduced by Smoluchowski J^l and extensively used in flocculation modelling |lfil 
1171 131 ITU 15], As also the bacterial growth © is compatible with the population balance 
structure, our model can be considered as an extension of PBM. 

In contrast to the bacterial growth dynamics the floe interactions (JHJ satisfy conser- 
vation of bacteria density, i.e., 

n —\ \ / floe interaction 



Attachment— detachment dynamics 

For a subset of the floe interactions (0, we are able to push the mathematical analysis 
further. In particular, we will restrict the interactions to the aggregation of a floe with 
a single bacterium, and the splitting off of a single bacterium from a floe. We call these 
processes attachment and detachment, respectively. These interactions can be written as 

F n + F\ — > F n+ \ with reaction rate <z„, 

F n+ i — > F n + Fi with reaction rate b n . (9) 
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The part of the dynamics corresponding to this type of floe interaction is 



dui 
~~dT 

du n 
~df 



= -2a\u\ - a n mu n + 2b\u 2 + &«-iu„ 

floe interaction n—2 n—3 

= a„_iUitt n _i - a„uiu n + b n u n+ i - 6„_itt n , n > 2. (10) 

floe interaction 



Since these equations are a special case of (jHJ), they satisfy again conservation of bacteria 
density. 



3 Fast flocculation dynamics 

The previous section introduced the model with © and (JHJ) or (|10(l . In this section 
we present an analytical model analysis, by assuming that the flocculation dynamics are 
much faster than the bacterial growth and the reactor dilution. This assumption will be 
relaxed in Sect.^] 



Separation of time scales 

Population balance models when used to describe flocculating bacteria, assume that the 
flocculation can be uncoupled from the other processes. It is argued that in the settling 
tank the substrate concentration s is sufficiently low to justify this assumption. The 
situation is however less clear in the reaction tank. Literature reports flocculation times 
of the order of 1 to 10 minutes |181 15]. to be compared with bacterial growth times, i.e. 
the inverse of the specific growth rate h(s), of 1 hour to 1 day and with retention times, 
i.e. the inverse of the dilution rate D, of a few hours to a few days. 

Even if the separation of time scales is not always satisfied in reality, it is interesting to 
investigate how our model behaves as this separation becomes infinitely large. Indeed, 
the model simplifies drastically and can be studied quite explicitly, as we will now show. 
Moreover, the reduced model can be considered as an approximation for the full model, 
as we will show in the next section. 

To make the separation in time scales explicit, we introduce a small parameter e > 0, 
du n ( du n \ 1 ( du n \ 

-TT = -TT I - D U„ + - — I . (11) 

Cir \ OC / bacterial growth 6 V QI / floe interaction 

Taking e — > 0, we introduce a sharp distinction between 

• the fast dynamics, consisting of the floe interaction, for times t ~ e, and 

• the slow dynamics, consisting of the bacterial growth and the dilution, for times 
t ~ 1. 



The idea now is as follows. On the short time scale, the system evolves to fast dynamics 
equilibria (u^f st ) which are parametrized by the total bacteria density x. On the large 
time scale, the system evolves on the manifold of these equilibrium distributions. As this 
manifold is one-dimensional and parametrized by x, we obtain autonomous dynamics for 
the biomass density x. 

To be more explicit, let us write down the dynamics for x by introducing © into Q). As 
the floe interactions conserve x, 

dx / du n \ ^ 

n —\ V / bacterial growth n— 1 
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Using © and JHJ), 

d °° 

— = ^ h n {s)u n - Dx. 

n—1 

At this point we use the separation of time scales. As we are looking at the slow dynamics, 
the distribution (u n ) will have reached its fast dynamics equilibrium (u^^x)) . Therefore, 

, 00 

= /i«(s)< ast (x)-Dx = h(s, x)x-D x, (12) 
n=i 

where we have introduced the specific growth rate h(s,x), 

00 

X>„(*R ast (*) 

h(s,x) = — . (13) 

X 

Note that h(s,x) depends both on the substrate concentration s and the biomass concen- 
tration x. The obtained specific growth rate is therefore density-dependent, as announced 
in the introduction. 



Uniqueness of fast dynamics equilibrium 

In the previous computation, we used the hypothesis that the fast dynamics reach the equi- 
librium distribution (u^ st (x)). This hypothesis can be justified by exploiting the analogy 
between molecules and chemical reactions on one hand, and floes F n and floe interactions 
(|7|) on the other. Indeed, the condition for chemical equilibrium can be translated to 

'■' '' for all in. n. 



with K m ^ n the equilibrium constant, independent of any density Uk- Obviously, not all 
these conditions are independent. For example, 

Fin -|- Fn ^ — F rn j rn 

Fm+1 g F m + F\ 

Fjn+1 ~t~ F n ^ — F. m + n +\ , 

and the equivalent computation for the equilibrium constants, 

jy- ^m+n **m+n+l " , m"l 7^ j> — 1 

-**-m+l.n — — **-m,n-**-Tn-\-n,l-ft-'m 1* 

To compute the fast dynamics equilibrium, it suffices to consider a basis of chemical reac- 
tions, i.e., a set of independent reactions from which the other reactions can be obtained 
by taking linear combinations. For our case, one such basis is given by the attachment- 
detachment interactions ©■ An even simpler basis is 

nF\ f± F n with equilibrium constant K n . 

The equilibrium conditions then read 

u n = K n u^, for all n > 2. (14) 

We now claim that for any total bacteria density x, there is only one distribution (u n ) 
satisfying (|14l) and Indeed, by l|14f> all u n with n > 2 are expressed in terms of u\. 
Introducing this in |JS}) 

oo oo 

x = nu n = nKnU™. 

n—1 n—1 

The right-hand side is a polynomial in u\ with positive coefficients, and vanishes for 
u\ — 0. As a consequence, for every positive x, there is a single u\ for which this equation 
is satisfied. Using (|14[) we obtain u n for all n > 2. This yields the unique distribution 

(Un St (x)). 
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Stability of fast dynamics equilibrium 

The previous analysis started from an equilibrium description for the floe interactions. 
The validity of such a description can be argued by noting that on the short time scale, 
our model reduces to a closed system of chemical reactions. Thermodynamics guarantees 
that such a system evolves towards a unique equilibrium, which can be computed in terms 
of the equilibrium constants K m ^ n I2()j . 

Nevertheless, it would be more satisfying to derive our conclusions directly from the kinet- 
ics JSJ for the floe interactions. However, this is only possible under additional assumptions 
on the floe reaction rates a rn . n and b mM . A sufficient condition is the so-called detailed 
balance condition, which essentially boils down to requiring the positivity of the entropy 
production I^Hj . We do not develop this argument further. 

It is interesting to note that in the case of attachment-detachment interactions 0, the 
uniqueness and the stability of the fast dynamics equilibrium can be shown explicitly. This 
is consistent with our previous remark, as the detailed balance condition can be shown 
to become empty for these interactions. In other words, the application of the separa- 
tion of time scales can be justified mathematically in the case of attachment-detachment 
interactions 0- 

First, we prove the uniqueness of the fast dynamics equilibrium for attachment-detach- 
ment. We have to find a distribution (u n ) such that the right-hand side of (|10|) vanishes. 
As these equations satisfy conservation of bacteria density, we can neglect one of the 
equations, say n = 1. Summing all the other equations, we obtain 



As this equation has the form of l|I4(l , the argument following the latter equation can now 
be repeated. 

Next, to prove the stability of the fast dynamics equilibrium (u^f st (cc)) for attachment- 
detachment, we make the additional assumption that the coefficients a n and b n have the 
same size dependence, i.e., 



a\U\Ui — b\u-2 = 



or 




This relationship allows us to express all the u n , n > 2 in terms of u\ , 




there exists a function n 



/„ such that a 



'71 



n 



and b n = bf ., 



Consider then the entropy-like functional 




Computing its time derivative under the fast dynamics 




n=l 



n=l 




oo 




Moreover, it vanishes if and only if 



for all n. 
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Together with the conservation of total bacteria density x, we conclude that H(t) is strictly 
increasing in time until it reaches its maximum at (u n ) — (u^f st (ir)). This establishes the 
announced stability property. 

Density-dependence of growth rate 

Under the assumption of fast flocculation, we derived the specific growth rate (|I3(1 . To- 
gether with 10 and {T1J , 

h(s,x) = with u^ at (x) = K n (uf st (x)) n . 

X 

For a = 1, we find h(s,x) = h(s). Indeed, the case a = 1 corresponds to the assumption 
that all the bacteria in the floes have the same substrate density s available. The specific 
growth rate is then substrate-dependent. On the other hand, when a < 1, we obtain a 
genuine density-dependent growth rate. 

The difference between substrate-dependence and density-dependence has important eco- 
logical consequences. Indeed, consider different bacterial species growing on a single sub- 
strate. A model with substrate-dependent growth rates predicts an equilibrium where 
only one species survives, whereas a model with density-dependent growth rates allows for 
an equilibrium where different species coexist. To establish the latter property, Ref. |?T] 
requires the map x i— > h{s,x) to be decreasing. We now present a proof of this property 
for our growth rate h(s,x). 

To simplify notation, we drop the superscript "fast" . As x <— > Ui(x) is increasing, it suffices 
to prove that h(s,x(ui)) > h{s,x(u'i)) if U\ < u' l7 or 

n ^ n 

nK n u™ ^ nK n u[ n 

n n 

Equivalently, 

^m a K m u^nK n u' 1 n -~J2m a K m u' 1 m nK n u^ > 0. 

m,n m,n 

As the terms m = n cancel out, this difference can be written as 

m<n 

Both factors in brackets are positive because a < 1, m < n and u\ < u^. This finishes the 
proof. Note also that the inequality is strict as soon as K m > 0, K n > 0, m ^ n. 



4 Slow flocculation dynamics 

The analysis in the previous section is based on the hypothesis that the parameter e in 
Eq. Ijlll) is small. This separation of the flocculation and the bacterial growth time scales 
seems to be well satisfied in the settling tank, but rather inadequate for modelling the 
reaction tank. In this section we investigate how our analysis can be extended without 
making this assumption. 

Recall the two dynamical systems we are considering. First, there are the full dynamics 
in terms of the densities u n of floes of size n. Together with the substrate dynamics, the 
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Table 1: Parameter values used in the simulations 



floe growth rate h n (s) 


h n (& 


) = h(s) n a with a = | 


bacterium growth rate h(s) 


h(s) 


0.2 s 

~ s+6 


attachment coefficients a n 


a n = 


- n a with a = | 


detachment coefficients b n 


K = 


: 0.1 n with a = g 


dilution rate D 


D = 


0.04 


inflow substrate concentration Sj n 


Sin = 


= 20 



system is given by 

dui , . 1 /dwi 

—— = -/ii(s)mi - Dili + - I -t— 

dt 6 V dt / fioc interaction 

du n If du n . 

h n -i(s)u n -i - h n (s)u n - Du n + - —t— n= 2,3, ... 



floe interacts 



ds °° 

^ = - / ] h n (s)u„ + D (si n - s). (15) 

n— 1 

Under the hypothesis of separate time scales, the full dynamical equations can be approx- 
imated by the reduced dynamics 

— = h(s, x)x — Dx 
ds 

— = -h(s, x)x + D(s- m -s). (16) 



We integrate both dynamical systems with parameter values given in Tab. ^ The floe 
growth rate behaves as h n (s) oc n a with a = |, corresponding to the surface-to- volume 
ratio of spherical floes. We consider attachment-detachment interactions (|10|) . with re- 
action rates a n oc n a and b n oc n a with a = |. The infinite sequence n = 2, 3, . . . of 
dynamical equations (|15|l is truncated at n = 300. The parameters of Tab. H and the 
initial conditions in the simulations are chosen such that this truncation yields a good 
approximation of the full dynamics. 

Fig. ^ compares the full dynamics l|15(l for different values of the parameter e with the 
reduced dynamics l)16[l. For small e, the solutions of l|15|) for different initial conditions 
converge rapidly (after a time of the order t ~ e) to each other. The solution of 116|l almost 
coincides with those of the full dynamics, indicating that the latter can be approximated 
as dynamics on the manifold of distributions (u l ^ st (x)). When e increases, the solutions 
for different initial conditions differ more and more. This indicates that there are no longer 
autonomous dynamics in the variable x, and thus no well-defined specific growth rate. 

We conclude that for larger values of e, the system cannot be described by a dynamical 
equation of the form l|16(l . Nevertheless, Fig. ^ shows that for all values of e, the different 
initial conditions lead to the same equilibrium. On the other hand, the equilibrium of 
(|16|l different from the wash-out solution, i.e. x ^ 0, satisfies h(s,x) = D. If we want 
the reduced dynamics to predict the correct equilibrium, the specific growth rate should 
satisfy this condition. In this way, we obtain a well-defined density-dependent growth rate, 
which we call the specific growth rate at equilibrium h° qul (s,x). 

Fig. plots the specific growth rate at equilibrium for different values of the parameter 
e. For small e, the specific growth rate at equilibrium coincides almost with the explicit 
formula JT3J. As e increases, the difference with JTSjl becomes substantial. Moreover, the 
function h cqnl (s,x) no longer factorizes as a function of s times a function of x. Note 
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Figure 1: Comparison between the full dynamics 1)15(1 for different e and the reduced 
dynamics 1(16(1. Parameter values of Tab. ^ were used. The full dynamics, shown in full 
line, were integrated for six initial conditions: three with x(0) = 20 (u n (0) = 20<5 ni i, 
w n(0) = 2<5 nj io and u n {0) = 0-2S n ,wo) and three with x(0) = 5 (u n (0) = 5S Ui i, u n (0) = 
0.5 5 nj io and u„(0) = 0.05 <5„.ioo). The reduced dynamics, shown in dashed line, were 
integrated with initial conditions x(0) = 20 and x(0) = 5. (a) e = 0.001; (b) e = 0.01; (c) 
e = 0.1; (d) e = 1. 
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Figure 2: The specific growth rate at equilibrium ft, equl (s,a;) as a function of the biomass 
density x for a fixed substrate density s = 6. Parameter values of Tab. ^ were used. The 
curves in full line correspond to, from bottom to top, e = 1, e = 0.1 and e = 0.01. The 
specific growth rate 1(13(1 is shown in dashed line. 
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Figure 3: Response of the full model I) 15)1 and the reduced model (|16|l to a step in the 
dilution rate D. (a) Excitation in the dilution rate D. (b) Reaction of the two models. 
Parameter values of Tab.n were used. The full line corresponds to the system 1)150 with 
e = 1. The dashed line corresponds to the system 1)16)1 with specific growth rate h eqnl (s, x). 
For both simulations, the initial condition was taken as the equilibrium for dilution rate 
D = 0.04. 



that all these curves are monotonically decreasing in the bacteria density x, as we proved 
explicitly for 1)131) . 

The reconstructed growth rates h oqm (s, x) can now be used to integrate (|16fl . i.e., 



d.r 

~d7 
ds 

dt 



= h eq ™(s ) x)x-Dx 
= -h cqui (s,x)x + D( Sil 



By construction, this model will tend to the same equilibrium as the full model 1)1 10 ■ To 
test how well it approximates the dynamics 1)11)1 , we perturb the system out of equilibrium 
and look at the resulting dynamics. As shown in Fig. ^ perturbations which disturb too 
heavily the floe size distribution cannot be correctly modelled by an equation of the form 
()160 . We therefore apply a perturbation in the dilution rate D, which acts similarly on 
the different floe densities u n . Fig. [3] shows that the reduced model predicts with rather 
good precision the reaction of the full system to this perturbation. 



5 Conclusion 



In this paper we investigated how flocculation influences the bacterial growth dynamics 
in a bioreactor. In the context of the activated sludge process, this coupling of physico- 
chemical and biological phenomena is mostly relevant for the reaction tank. In particular, 
we studied the possibility of an effective model on the level of the biomass density, without 
explicitly taking flocculation into account. 

Such an effective description is only possible when the flocculation dynamics are sufficiently 
fast compared to the other processes. In this case, the specific growth rate, which for 
isolated bacteria depends only on the substrate density, gains an additional dependence 
on the biomass density. It is interesting to note that such a density-dependent growth rate 
has recently been proposed as a mechanism to explain the coexistence of many bacterial 
species growing on a limited number of substrates. We will investigate the link between 
flocculation and species coexistence in a forthcoming contribution. 

When the flocculation dynamics have time scales comparable to the bacterial growth, 
the details of the floe size distribution do affect the global system dynamics. In that 
case, dynamics autonomous in the biomass density do not exist, and the notion of specific 
growth rate is ill-defined. However, if the reactor evolves such that the floe size distribution 
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remains equilibrated, it makes sense to define a specific growth rate at equilibrium. We 
showed in a simple example, that such a growth rate, which is again density-dependent, 
can yield an accurate description of the system dynamics. 
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